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Abstract 

The quest to build a mass formula which have in it the most relevant microscopic contributions 
is analyzed. Inspired in the successful Duflo-Zuker mass description, the challenges to describe the 
shell closures in a more transparent but equally powerful formalism are discussed. 
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I. INTRODUCTION 



Understanding nuclear masses provides a test of our basic knowledge of nuclear structure. 
Its accurate knowledge is relevant for the description of various nuclear and astrophysical 
processes |l|], Though great progress has been made in the challenging task of measuring 
the mass of short-lived nuclei which are far from the region of stable, naturally occurring 
isotopes, theory is needed to predict their properties and guide experiments that search, for 
example, for regions of increased stability [2]. 

Advances in the calculation of atomic masses have been hampered by the absence of a 
full theory of the nuclear interaction and by the difficulties inherent to quantum many-body 
calculations. There has been much work in developing mass formulas with both microscopic 
and macroscopic input, on one side, and on the derivation of masses in a fully microscopic 
framework, on the other [3J]. The most successful approaches are the Finite Range Droplet 
Model (FRDM) the Skvrme and Gogny Hartee Fock Bogolyubov (HFB) [5], and the 
Duflo-Zuker mass formula p-lSj. They allow for the calculation of masses, charge radii, 
deformations, and in some cases also fission barriers. They all contain a macroscopic sector 
which resembles the Liquid Drop Mass (LDM) formula, including volume and surface terms, 
the Coulomb interaction between protons, Wigner and symmetry terms, linear and quadratic 
in the neutron excess N — Z, and pairing. They also include deformation effects in different 
ways. HFB calculations have succeeded in going through the 1 MeV barrier, which until 
very recently seemed unsurmountable, and have achieved RMS deviations smaller than .6 
MeV 9], competitive with the published fits obtained with the FRDM, which are also being 
improved. On the other hand, the astonishing RMS around 0.3 MeV obtained by the DZ 
model 15 years ago sets the standard and provides the most robust predictions [3! \l0\. 

The Liquid Drop Mass (LDM) formula captures the macroscopic features of the mass 
dependence on the number of neutrons N, of protons Z, and on its mass numbers A = N+Z. 
Shell effects refer to the differences between the experimental binding energies Uj and the 
^DM predictions. Based in the remarkable success of the Duflo-Zuker (DZ) mass formula 

in n 

6-8], the nuclear monopole Hamiltonian [12] allowed the description of the energy spectra 
of particle and hole states and particle- hole gaps on double magic cores jl3 |. Alternatively, 
there have been various attempts to describe the shell effects through a simple inclusion 
of linear and quadratic functions of the number of valence nucleons. They improve the 



description reducing the root mean square deviation (RMS) by half 

mm When 2-, 3- 

and 4-body terms are included, following the simplest DZ approach, the shell corrections for 
nuclear masses and radii can be calculated fl6 |. 

The DZ mass model provides an attractive combination of simplicity and microscopic 
components. Since its initial formulation jcj-8|, there have been efforts to communicate its 
philosophy [171. [lSfl. A detailed description of the simpler DZ10 model, which contains the 



basic ingredients, is available 



19] . In the present contribution we present some general ideas 



for building a microscopic mass formula and speculate about the way in which they could 
be implemented. 

II. THE MACROSCOPIC SECTOR 



The macroscopic liquid drop model employed in the present analysis is a slightly refined 
version of the Liquid Drop Model. It contains a volume term, proportional to A, and a 
surface term, proportional to A 2//3 , a Coulomb term 

Z(Z -1) + 0.76[Z(Z -l) 2 / 3 ] 
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which includes the charge radius 
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instead of the mass radius 
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They both include the dependence of the nuclear radii with the isospin T = \N — Z\/2. 
The asymmetry term depends on the total isospin T 

4T(T+ 1) 
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The surface asymmetry is 
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The Pairing term includes corrections of order ^ and scales as l/r ~ 1/A 1 / 3 . It has 
different functional forms depending on whether N and Z are even or odd, listed in Table [B 
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TABLE I: The different expressions employed for the pairing contribution ep a %r are listed in the 
fourth column. The six cases are classified with the parity of N and Z, listed in the first two 
columns, and with N larger or smaller than Z, third column. 

With this six terms the macroscopic binding energy is built 

BE LDM = ai A - a 2 A 2/3 + a 3 e Cou i 

— 0-4 S-Asym + «5 esAsym + 0-6 ^Pair- (4) 

The coefficients a, have units of MeV. They were selected to minimize the root mean 
square deviation (RMS) when the predicted binding energies BE th (N, Z) are compared 
with the experimental ones BE exp (N, Z), reported in AME03, modified so as to include more 
realistically the electron binding energies as explained in Appendix A of Lunney, Pearson 



and Thibault 
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RMS _ I E [BE cxp (N, Z) - BE th (N, Z)f \ ^ (g) 
N nuc i is the number of nuclei for N, Z > 8. The minimization procedure uses the routine 



Minuit 
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The shell effects, i.e. the differences between the binding energies evaluated with the 
macroscopic LDM, Eq, (j4]), and the measured ones taken form AME03, are plotted in Fig, 
[1] as functions of N, Z, A and N — Z. They display information for the 2149 nuclei with 
N,Z > 8. To visualize the binding energies when only protons or neutron shells are open, 
we present also results obtained employing only the masses of the 185 nuclei in AME03 
which are semimagic, i.e. those with = 14, 28, 50, 82, 126 or Z = 14, 28, 50, 82. The extra 
binding associated with the shell closures are clearly seen. These effects are the target to 
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FIG. 1: Shell effects, in MeV, for the Liquid Drop Model, for the whole set of nuclei in AME03. 
and for the semimagic nuclei (heavy black dots). 

be described by the microscopic terms. Notice that the closed shell peaks are connected 
by smooth parabolas, evident in the plot for the semimagic, while the bottom regions are 
flat, a deformation effect taken explicitly into account in the DZ construction jf^O]. We can 



also observe that there is very little dependence on the isospin. The fitted coefficients are 
presented in the second column of Table [Til The RMS of the fit is 2.41 MeV. 



III. THE SHELL CLOSURES 

One special aspect of the DZ approach is its ability to build the volume and surface terms 
starting from the monopole interaction. The basic building blocks of the master terms are 
weighted sums over the occupied orbitals 

PvF PttF 
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where the sums are defined for neutrons {y) and protons (tt), and they run over all occupied 
p u , p w orbitals up to the Fermi level p u p, p v p with occupations n v , n % and degeneracies D u , D w . 
They are employed to construct the master term: 

MA = ^-(e ln +e lv )\ (7) 
This expression represents the dominant monopole component of the two-body interac- 



tion. The weights -7= were fitted from the dominant term of realistic interactions 17H19] . 

Performing the sums for closed harmonic oscillator shells p, with degeneracies D p = 
(p + l)(p + 2), it is not difficult to obtain the scaling of MA, which is proportional to the 
mass number A. In this way 
Its precise asymptotic behavior is 

MA a = ^(N^ + Z^) 2 . (8) 

The difficulty with this approach is that the shell closures are the ones associated with 
the 3D harmonic oscillator (HO), and the mass model has to include other contributions 
which displace the closures form the HO to the extruder-intruder (EI) scheme shown in Fig. 

12 

The strategy which we want to explore here is based in the fact that expressions ([6]) 
and ((71) can be redefined, with the index p u ,p n now referring to the EI major shell with 
degeneracies D U1T = D Vv n + 2. The fits of the monopole part of the two-body interaction 



cannot distinguish between the above mentioned one and this proposal 20| . While it can be 



argued that the harmonic oscillator basis is the simplest and natural one to be employed in 
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FIG. 2: Harmonic oscillator shells and extruder - intruder (EI) subshells. 

the microscopic description of the atomic nuclei, in what follows we show that the present 
proposal introduce in a very simple way the shell closures of the Extruder-Intruder single 
particle scheme, at iV = 14, 28, 50, 82, 126 and at Z = 14, 28, 50, 82. 

We have explored two different ways to include the shell effects. In the first one we add 
one term to the binding energy, proportional to MA — MA a . Its inclusion allows for a fit 
with an RMS of 1.46 MeV, a significant reduction for a single term. The fitted coefficients 
are listed in the third column of Table [TTJ 

Fig. |3]show MA — MA a as function of N, Z, A and N-Z for the whole set of nuclei included 
in AME03 and for the semimagic nuclei. The peaks associated with the shell closures are 
clearly seen, but the parabolas are now partially truncated. In between these peaks, new 
parabolas have emerged, associated with open-shell nuclei. 

Fig. H] displays the shell effects, i.e. the differences between the experimental binding 
energies and those obtained with the fit. The shell effects are still there, but diminished, 
with maximum differences around 4 MeV instead of 12 MeV for the pure LDM. 

The alternative approach is closer to DZ: to replace the volume and surface terms by MA 
and MA/r. To get a glimpse of its significance, in Fig. [5] we have plotted the differences 
between the binding energies obtained employing the master terms and the LDM. It is clear 
that the master terms are building the shell effects. The image for the semimagic nuclei is 
particularly explicit. 

In Fig. [6] the shell effects for the LDM with the volume and surface terms replaced by 
the master terms are presented. The master terms are somehow overestimating the shell 
closures, which now appear with negative sign. The parabolas are there, and the RMS of 
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TABLE II: The fitted coefficients, in MeV, with their associated mean and RMS errors. 
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the fit is 3.33 MeV, with the coefficients listed in the fifth column of Table HT1 

Up to here we have explored the inclusion of the EI shell closures trough the use of the 
master terms with EI major shells. The next step would be to include two, three and four 
body terms, of the type employed in DZ. All of them, being shell corrections, should scale as 
r w A 1 / 3 , and vanish at shell closures. It means that only valence particles contribute 
to them. At variance from the master terms, there are not sums here. They can be built 
with operators which vanish at closed shells, i.e. with zero valence particles or valence holes. 
They have, for both protons and neutrons, the generic form 

The details can be found in Ref. [l7|, 1 191 ] . 

The other challenge is the inclusion of deformation effects. This was successfully per- 
formed in DZ exploring different occupation numbers in the intruder and normal parity 
subshells. Inspired in the Nilsson diagrams, deformed nuclei were assumed to have four pro- 
tons and four neutrons removed form normal parity orbitals and occupying the next intruder 
orbital. From this point of view, the DZ model is a functional of the occupancies. One can 
visualize different models with increased complexity with DZ on one side and HFB on the 
other. 
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FIG. 4: Shell effects, in MeV, for LDM + MA - MA a (RMS=1.46MeV), for the whole set of nuclei 
in AME03 and for the semimagic nuclei (heavy black dots). 
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IV. CONCLUSIONS 



In the present contribution we have explored the possibility of describing the shell closures 
employing master terms of the type proposed by Duflo and Zuker, but with major shells 
built in the EI scheme. It seems more promising to add a simple master term with its 
asymptotic behavior removed that to employ the full master terms instead of the volume 
and surface terms of the LDM. If this approach, when combined with the inclusion of two-, 
three- and four-body terms and allowing for different occupation numbers, can lead to an 
useful mass formula is indeed an open problem. 
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FIG. 5: Differences between the binding energies obtained employing the master terms and the 
LDM, for the whole set of nuclei in AME03 and for the semimagic nuclei (heavy black dots), in 
MeV. 
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FIG. 6: Shell effects, in MeV, for LDM with master terms, for the whole set of nuclei in AME03 
and for the semimagic nuclei (heavy black dots). 
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